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I I Abstract 

rS^ This paper has evolved out of our previous work on static stress trans- 

Ph fer, where we used the full-space elastostatic Green's tensor to compute 

f-\ the Coulomb stress transfer impact of the Landers earthquake on the Hec- 

flj tor Mine event. In this work, we use the elastostatic Green's tensor for an 

(JJQ arbitrary layered Earth model with free-surface boundary conditions to 

ry\ study the impact of elastic heterogeneity as well as source-fault slip and 

^ geometry on the stress transfer mechanism. Slip distribution and fault 

* ^ geometry of the source have a significant impact on the stress transfer, 

^►s. especially in case of spatially extended triggered events. Maximization 

ri of the Coulomb stress transfer function for known aftershocks provides a 

O I mechanism for inverting for the source event slip. Heterogeneity of the 

elastic earth parameters is shown to have a sizeable, but lower-magnitude, 

^^ impact on the static stress transfer in 3D. The analysis is applied to Lan- 

►^ ders/Hector Mine and 100 small "aftershocks" of the Landers event. A 

I computational toolkit is provided for the study of static stress transfer for 

rvj arbitrary source and receiver faults in layered Earth. 

00 



1 Introduction 

^^ In this work we analyze the sensitivity of the static stress transfer mech- 

, anism (|6]) to the heterogeneity of the elastic earth model, source fault 

. . geometry and slip. Our approach to computing the change in the Coulomb 

Failure Function (CFF - see |15p due to a triggered stress change is sim- 
ilar to e.g. [6] and [12], and the one used in [8]. The results are ap- 
, ■ plied to analyzing a possible Coulomb stress transfer impact of the 1992 

^ M=7.3 Landers earthquake on the M=7.1 1999 Hector Mine event and 

100 smaller regional "aftershock" events. A connection between Landers 
and Hector Mine was conjectured soon after the occurrence of the second 
event; a static-stress impact causality analysis was presented in !8j and 
will be compared to our results. 

We numerically compute the elastostatic Green's tensor for an arbi- 
trary layered Earth model, and use the obtained tensor in the double- 
couple representation of displacement due to slip on a finite fault (|lljl. 



As in our previous work, the obtained displacements are used to compute 
the induced change in the stress field and OFF. To validate our approach, 
first we test our heterogeneous displacement modeling algorithm on a se- 
ries of source events with varying geometry and magnitude. 

In the application to Landers and Hector Mine, we use three differently 
oriented receiver fault planes and a multi-fault source model. Maps of 
CFF, normal and shear stress change at all points of the receiver faults 
are generated, and the impact of input parameters on the receiver slip 
likelihood is studied. 

Results of both stress transfer modeling and slip inversion algorithms 
are sensitive to the accuracy of input data for the source and receiver 
events (e.g. fault geometry). The focal mechanism inversion data for 
the Landers and Hector Mine events, provided by the Southern Califor- 
nia Earthquake Data Center (see [9]), fails a simple sense check: e.g., 
the data indicates that the Landers "fault" has the down-dip azimuth of 
165° - this disagrees with the known geometric features of the Johnson 
Valley, Landers, Homestead Valley, Emerson and Camp Rock faults that 
experienced temporally and spatially distributed slips during the Landers 
earthquake (see [IS])- The reason for such discrepancies is that the focal 
mechanism inversion "reconstructs" each event as a slip on a single fault 
based on early arrivals. More accurate fault data can be obtained from 
alternative inversion/observation methods that take into account regional 
displacement measurements. In this study, we use the fault orientation 
and earthquake nucleation data for the Landers event presented in |13j . 
Likewise for the Hector Mine event, rather than using the SCED-provided 
single-fault inversion parameters, we use the parameters estimated in [5], 
[3] and [10]. In particular, we compute the Coulomb stress transfer for 
three fault orientations/locations associated with the nucleation of the 
Hector Mine event (see [5]). The 100 small Landers "aftershock" events 
used in the analysis are taken from the SCED small-event database ([2]). 
Along with using the original fault parameters from the database, we al- 
low for the uncertainty in the parameters and create three "aftershock" 
datasets made up of events local to Landers, spread around the entire rup- 
ture, and a decimated version of the latter. The direct and inverse stress 
transfer analysis is performed for all three datasets and Hector Mine. 

The entire heterogeneous displacement modeling/ stress transfer frame- 
work has been implemented using MATLAB for easy reproducibility. A de- 
scription and brief user manual for the provided tools are included in 
Appendix A. 



2 Method 

We model the source earthquake as a series of spatially- distributed slips 
on finite faults in 3D. In the modeling part of this exercise, the source 
events are specified by 3D fault geometry (coordinates of the fault center, 
fault length, width, strike, dip and rake) and moment magnitude M^ 
for each fault. Once the elastostatic Green's tensor is obtained, the 3D 
displacement field due to the slip on the source fault is modeled using the 



VoHerra equation (see e.g. [11]): 



Mfc(x) = j ^is 
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%fe(x,y) 
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Tlj CliJy . 
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In (fTl), we have omitted the tensile (diagonal) components of the tensor 
in the square brackets as we assume that the slip is parallel to the fault 
surface, the integration is carried out over all faults UF, and (7fc(x,y) is 
elastostatic Green's tensor - i.e., the displacement along axis k at point 
X due to a unit force applied along axis i at point y (see e.g. |11)1. 
Variables p, Si, rij have the usual meaning of the shear modulus, fault slip 
and external normal vector components. Note that the tensile components 
of stress can be easily included in ([l| if required for e.g. modeling tensile 
cracks. 

For a homogeneous half-space earth model with the free-surface bound- 
ary conditions, the elastostatic tensor was given by Mindlin in /7 . In case 
of a fixed and known beforehand number of elastic layers, the Green's 
function can be found analytically using e.g. the method of images (|11|. 
[4]), or the analytic form of the propagation/stiffness matrix method ([4]). 
If the number of layers are unknown, the numerical form of the propa- 
gator matrix may be used (e.g. [I]). However, in this paper we use a 
different approach based on numerically solving the system of elastostatic 
equations with the free-surface boundary conditions &t z — and assum- 
ing zero co-seismic displacement at some maximum depth (20km in our 
tests). More specifically, we solve the elastostatic wave equation in the 
Navier form 

^("^)^" + r^M^ a^ a^ ^ °' ^ = ^' 2' '• (2) 

After Fourier-transforming (pi with respect to the two horizontal variables 
we obtain 
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where A = \{x-j,) = 2/i!^/(l — 2v) is the Lame coefficient. Rather than 
solving PI for both displacements and stresses, we discretize in z and 
complement the resulting system with the condition of continuity of nor- 
mal tractions ([11]) between discretization intervals. Once the elastostatic 
tensor is computed, the regional stress change is computed from forward- 
modeled displacements ([l| as 



Ao-fc = /i 



dui 
dxk 



+ 



duk 
dxi 



L-2!/ 



divu5fe. 



(4) 



where tensile stress is positive. Coulomb stress transfer (ACFF) at the 
receiver fault is then computed form Q, taking into account the stress 
sign conventions, as follows (see [IS]*): 



A(T„ — n,. AaUr 
At — Sr Aarir. 



ACFF = Ar + /i/ (Acr„ - APp) = Ar + M/Acr„ (1 - B) , (5) 



where n^, s,- are the external unit normal and unit slip (rake) vectors for 
the receiver fault, Aa„ is the normal stress change (positive in tension), 
and At is the shear stress changaj induced by the source events. The 
Skempton coefficient B relates the change in the confining stress to the 
resulting change in pore pressure f[14p and is a source of uncertainty in our 
analysis. In this work we assume B ~ 0.5 based on the fact the range of 
values for solid rocks is ~0.5-0.8 ([E]). The coefficient of fault friction fif 
typically ranges between 0.5 — 0.9 and is depends on gouging/brecciation 
within the receiver fault. We attempt to estimate this parameter based 
on plausible stress transfer scenarios for the analyzed source and receiver 
faults. 



Equations (Tl5l provide a mechanism for mapping known slips on the 



source fault to the stress change on the receiver fault, and form the basis 
of the static stress triggering analysis. We additionally suggest to solve 
an inverse problem for identifying the slips on the source fault from the 
receiver fault parameters. In our proposed inversion method we do not 
use values of the stress change at the receiver faults but rather attempt 
to estimate the source slip from the assumption that the user-specified 
spatially-distributed receiver faults have been encouraged to slip by the 
source events. More specifically, if A,B,C are linear operators that map 
source slips to Acr„, At, ACFF 



Aa„ = As 
At = Bs 
ACFF = Cs, (6) 



then we formulate the following linear programming problem for comput- 



^in the direction of the receiver slip (rake) 



all grid points of all receiver faults 

[11 . . . 11] Cs — ^ max, 

grid of receiver fault # k 

[0...0 11... TT 0...0]As > 0, 

grid of receiver fault # k 

[0...0 11... TT 0...0]Bs > 0, 

grid of receiver fault # k 

[0...0 11... TT o...o]Cs > 0, 

li{Sinj + Sjni)dS = M^j, 

s > 0, (7) 

where k ranges across all receiver fault indices, / ranges across all source 
faults _F1 , s is a vector made up of non-negative slip along respective rake 
direction at each grid point of each source fault, M^j is the moment ten- 
sor associated with the source event (slipping fault) No I. Equations pi 
mean that we maximize the total change of the Coulomb Friction Func- 
tion ^ ACFF for all receiver faults under the constraints that the to- 
tal Acr„, At, ACFF for each receiver faulPj'is positive. Additionally, the 
slips are constrained by the known moment tensor and the obvious re- 
quirement that the slips are non-negative in the rake direction. The three 
stress-related constrains can actually be reduced, depending on whichever 
slip-encouragement mechanism is likely relevant (e.g., shear advancement 
or un-clamping) . This is especially useful with uncertain friction and 
Skempton coefficients. 



3 Implementation 



The elastostatic tensor is computed on a 3D grid for a number of source 
depths^ In all of our experiments we use the geological coordinate system 
with the axis x pointing North, axis y pointing East and z pointing down 
(depth). Negative vertical displacements correspond to uplift while posi- 
tive ones, to subsidence. The Green's tensor is computed for two modeling 
cases: 

1. 128 X 128 X 20 km on a 64 X 674 x 20 grid, with source depths (index 
of 2) ranging from 2 to 20; 

2. 256 X 80 X 20 km on a 256 x 80 x 20 grid, with source depths ranging 
from 2 to 20. 

The delta function was modeled as a Gaussian with the standard deviation 
of 3 and 1.5 km, respectively. The second grid has been chosen so as to 



■^integrated over the fault surface 

•^The Green's tensor for a layered Earth is invariant with respect to planar translations 
of both source and receiver hence the source can always be assumed to have x = 0,y = 0; 
however, source and receiver depths must be independent. 



encompass just the Landers/Hector Mine geography (see e.g. [13]). Local 
fault grids have x pointing in the strike direction, y down dip and z into 
the foot walj_] 

In our analysis of sensitivity to elastic heterogeneity, we use two earth 
models: homogeneous /x = 31 GPa and a horizontally layered one with 
more compliant top layers but a similar average shear modulus over seis- 
mogenic depths of ~15 km (see figfl}. Poisson ratio is assumed constant 
in both cases, equal to 0.2. The shear modulus and (computed) Lame 
constant are smoothed with a 4-km rectangular filter before solving S - 
see fig IT] 

The elastostatic Green's tensor computation is performed by the 
elastostatic_green.m MATLAB module, described in Appendix A. 





Figure 1: Homogeneous and layered Earth model. Smoothed model. 



A component of the Green's tensor computed on the square grid is 
shown on fig [2] Note that the more compliant top layers resulted in a 
smaller uplift in the compressive area to the north of the source. Figures 
[3j[4j[5j[6]present the results of forward-modeling surface displacements due 
to the Landers earthquake using the elastostatic Green's tensor computed 
on the 256 x 80 x 20 grid. The Landers is modeled as three slips along 

• the combined Johnson Valley/Landers fault (centered at 34.225°N, 
116.425°W, depth 9 km, strike -5°, length 30 km, width 14 km); 

• Hempsted Valley fault (centered at 34.425°N, 116.475°W, depth 9 
km, strike -26°, length 25 km, width 14 km); 

• Emerson/Camp Rock fauh (centered at 34.575°N, 116.575° W, depth 
9 km, strike -40°, length 30 km, width 14 km), 

with the constant dip of 90° and rake of 180°. The total Landers moment 
is distributed as 18%, 55%, 27% among the three "finite- plane" events ac- 
cording to the estimated spatial distribution of the slip over fault surfaces 
(see Fig 5 of [13]). 

Note that since the event is modeled as a strike slip close (2 km) but 
below the surface, the resulting behavior is consistent with the presence 
of more compliant top layers (see figp][6|: a smaller shear modulus on 
the fault near the top means greater slips that translate (across the thin 
2 km top layer) to greater horizontal displacements on the surface fig p] 
And a greater compliance of the medium near the surface results in a 



*e.g. for a strike-slip fault striking Nortii tlie local z axis in the fault coordinate system 
points West 



smaller uplift/subsidence (fig El. Note the uplift in the compressional 
areas and subsidence in the tensile regions (fig SlSl. Additional exper- 
iments for simple strike-slip and thrust faults indicate good qualitative 
agreement of the results of our modeling algorithm with the expected 
medium behavior ([11]). Surface displacement modeling is performed by 
the model_surface_displacements.m MATLAB module, described in Ap- 
pendix A. Note that this module is not required for the stress transfer 
computations and is provided simply for validating the elastostatic ten- 





Figure 2: Green's Tensor Component: vertical displacement due to unit force 
in Northern direction, E 90° strike, homogeneous vs layered Earth. 



Vertical Surface Displacements 



Vertical Surtace Displacements 
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Figure 3: Landers Forward Model, vertical displacement, homogeneous vs lay- 
ered Earth. 

Both forward and inverse stress transfer problems ([5| and m are 
solved by the transfer .stress. m module. 



4 Stress Transfer 



Fig Is] shows a regional stress map with some pre-Landers events indicated 
by dots. Fault and nucleation data for the Landers event ([13]) and Hector 
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Figure 4: Landers Forward Model, vertical displacement contour, homogeneous 
vs layered Earth. 
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Figure 5: Landers Forward Model, horizontal displacement directions, homoge- 
neous vs layered Earth. 
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Figure 6: Landers Forward Model, horizontal displacement, homogeneous vs 
layered Earth. 



Mine ([5]) place the source and receiver faults in block 72 of the map, with 
the Hector Mine faults situated to the North East of Landers. 

Stress tensor inversion studies for this region indicate the principal 
maximum regional stress directions near the block 72 ranging from -10° 
to 5° measured from North. The indicated stress directions are consistent 
with the Hector Mine faults being in a state of failure equilibrium (see 
[15|). as the estimated strike directions of the Hector Mine fault (-35° and 
-15° - see [5]) position those faults within « 25 — 30° of the maximum 
stress axis. Assuming an approximately strike-slip faulting regime, this 
range of angles corresponds to values of the friction coefficient between 
0.6 ^ l/tan2 x 30° and 0.8 ^ l/tan2 x 25°. 



Landers and potential aftershocks 
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Figure 7: Possible Landers aftershocks used in the inversion: (a) clustered 
around Landers (b) more uniformly distributed around the entire rupture. Cen- 
ters of Landers/Johnson Valley, Homestead Valley, Emerson/Camprock seg- 
ments of the rupture marked. 

We have randomly selected 100 events from the SCED database ([2 ) 
within 18 km from the center of the combined Landers/ Johnson Valley 
faults ('local dataset') and spread uniformly around the entire Landers 
rupture ('regional dataset') - see figlT] Additionally, we solve the direct 
and inverse stress transfer problems for a decimated regional dataset con- 
sisting of 20 randomly picked events. The results of stress transfer compu- 
tation for the local dataset are shown on fig |10|ll|12| The corresponding 
plots for the complete and decimated regional datasets are on fig |16|17|18| 
and fig |13|14|15[ respectively. In all cases majority of faults appear to be 
encouraged to slip through the advancement in shear and/or un-clamping. 
Green circles indicate that the corresponding quantity advances the slip, 
while red crosses denote the oppositcrl Due to the proximity of these 
aftershocks to Landers, elastic heterogeneity (e.g. compliance of the top 
layer) does not have a significant qualitative effect on stress transfer. How- 
ever, elastic heterogeneity becomes important at greater distances when 
the stress change is <^ 10 Bar. We will demonstrate the effect of hetero- 
geneity by computing the stress transfer on the receiver faults for Hector 
Mine. 



^symbol sizes are proportional to the corresponding values 




Figure 8: Regional Fault Map. 



The Hector Mine is included in a separate dataset and modeled, fol- 
lowing [S|, as three faults shown on fig lo] from [S', using the following 
parameters: 
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where x, y, z are the coordinates of respective fault centers with center of 
the Landers/Johnson Valley fault at the origin. Note the nucleation into 
faults 1 and 2 (see fig|9| that appear to spatially overlap if viewed e.g. in 
the eastern direction. 

Results of the stress transfer calculations indicate that un-clamping 
appears to be the dominant slip-encouraging effect of Landers on Hector 
Mine (see fig 19|20|23|24[ | . Stress transfer is computed on local fault grids 
of 7 X 8, 10 X 8 and 7x8 points. Note that due to the spatial overlap 
of faults 1 and 2 (see fig ^ and relative sizes of the two faults (see the 
parameters above) the northern half of fault 2 (i.e., the right half along 
the strike axis) can be expected to exhibit roughly the same properties 
as the entire fault 1. As can be seen from e.g. fig |19|21|23| this is indeed 
the case, with the difference simply due to the different strikes (15° and 
24°). The shear stress appears to drop throughout the Hector Mine rup- 
ture (fig 21|22 I, while a significant un-clamping is predicted on faults 1 
and 2 (fig 23|24|. Coulomb transfer appears negative almost everywhere 
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Figure 9: Hector Mine faults with the previously estimated slips, from [5]. 



except at the "junction" of faults 1 and 2 - i.e., where the maximum slip 
has been predicted along the Hector Mine rupture (see fig [9] and [5]). It is 
plausible that Landers statically triggered Hector Mine by stress transfer 
to this point that happens to coincide with the hypocenter ([5], [8]). A 
slightly higher friction 0.8 results in a greater slip encouragement, how- 
ever, does not alter the qualitative picture. Note that our results are in a 
good qualitative agreement with ^ (compare with plate 2 of [8]) and our 
previous estimates in HW7. 

The results indicate a significant sensitivity to the elastic heterogeneity 
of the medium. Stress plots appear to be shifted down for the heteroge- 
neous test, consistent with equal average moduli but more compliant top 
layers in the heterogeneous case. Note the resulting quantitative impact 
e.g. normal stress drop on fig |23|24[ This indicates the importance of 
accounting for elastic heterogeneity in stress transfer analysis. 



5 Slip Inversion 



As was demonstrated in the previous section, the estimated stress transfer 
is not always such as to encourage slip, hence some constraints in ([7| will 
be violated. Including all constraints may render the problem inconsis- 
tent and the solution infeasible for the optimization algorithm (Simplex 
implemented in MATLAB linprog). Therefore we conduct two experiments: 
invert using all the constraints, or use only consistent constraints obtained 
from direct modeling, allowing us to assess the impact of dataset size on 
the accuracy. 

Fig 25|26 demonstrate inversion results for the (three idealized) Lan- 
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Figure 10: Coulomb transfer for the local dataset of Landers aftershocks (20 
events), homogeneous vs layered Earth. 
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Figure 11: Advancement in shear for the local dataset of Landers aftershocks, 
homogeneous vs layered Earth. 



12 



E "^ 





-12 -10 -8 -G -4-2 2 



-12 -10 -8 -6 -4-2 2 



Figure 12: Unclamping (green) for the local dataset of Landers aftershocks, 
homogeneous vs layered Earth. 
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Figure 13: Coulomb transfer for the decimated regional dataset of Landers 
aftershocks (20 events), homogeneous vs layered Earth. 
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Figure 14: Advancement in shear for the decimated regional dataset ol Landers 
aftershocks, homogeneous vs layered Earth. 
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Figure 15: Unclamping (green) for the decimated regional dataset of Landers 
aftershocks, homogeneous vs layered Earth. 
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Figure 16: Coulomb transfer for Landers aftershocks, homogeneous vs layered 
Earth. 
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Figure 17: Advancement in sliear for Landers aftershocks, homogeneous vs lay- 
ered Earth. 
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Figure 18: Unclamping (green) for Landers aftershocks, homogeneous vs layered 
Earth. 
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Figure 19: Coulomb transfer, Hector Mine faults 3,2,1, /i — 0.7, homogeneous 
vs layered Earth. 
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Figure 20: Coulomb transfer, Hector Mine faults 3,2,1, /i = 0.8, homogeneous 
vs layered Earth. 
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Figure 21: Shear Stress Transfer, Hector Mine faults 3,2,1, fi — 0.7, homoge- 
neous vs layered Earth. 



17 





M ill 



i 






I 



i 



Figure 22: Shear Stress Transfer, Hector Mine faults 3,2,1, /i — 0.8, homoge- 
neous vs layered Earth. 
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Figure 23: (Un) clamping. Hector Mine faults 3,2,1, /x = 0.7, homogeneous vs 
layered Earth. 
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Figure 24: (Un)clamping, Hector Mine faults 3,2,1, /i = 0.8, homogeneous vs 
layered Earth. 



ders faults using the local dataset (i.e., events clustered around Landers). 
Despite the events being clustered around source fault 1, inversion indi- 
cates the most significant slips on source fault 2 (i.e.. Homestead Valley) 
consistent with our earlier estimate that this part of the rupture accounts 
for more that half of the total Landers moment. This is confirmed by the 
inversion of the decimated regional dataset consisting of only 20 faults: 
fig |27|28] again indicate the greatest moment density on the source fault 
2. Including the entire regional dataset produces qualitatively similar re- 
sults of fig |29|30] Of course, this is an obvious consequence of the equality 
constraint in (JTl) for the total moment on a source fault. However, an- 
other, more subtle, feature of all these plots is that fault 2 is predicted to 
have small slips roughly in the interior of the southern half of the fault, 
with the slips increasing between faults 1 and 2. Taking advantage of the 
spatial distribution of Hector Mine, we may expect inversion of Hector 
Mine to produce a "more accurate" prediction of slips on fault 2. Fig |31| 
demonstrate the result of solving for Landers slips using only un-clamping 
constraints for Hector Mine. The results appears to be in a qualitative 
agreement with the estimated Landers slips (see Figure 5 of [13]), cor- 
rectly identifying the low-slip area 7-10 km from the southern tip of the 
Homestead Valley fault, starting at a depth of ~7 km. While this result 
should not be over-interpreted, it may indicate inverse stress transfer as 
a viable mechanism for slip inversion. Unlike the traditional methods of 
inverting subsurface slips from the observed surface displacements (e.g., 
GPS or INSAR, [11]), inverse problem nl estimates the slips from the 
assumption that the source fault triggered a given set of events. Such in- 
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version is inevitably qualitative by nature but still may prove useful as 
e.g. an additional constraint or regularization for a more traditional in- 
version method. Further study of inverse stress transfer in applications to 
large-scale triggering problems may help assess its usefulness. For exam- 
ple, fig |32] shows the result of inverting the slips on Landers fault 2 using 
the decimated regional dataset with all constraints, on a high-resolution 
fault grid. The result is biased by the inconsistent constraints but appears 
not unreasonable in the context of known Landers slips ([IS]). A simi- 
lar experiment could not be carried out for Hector Mine as the receiver 
because of the numerical cost of solving this problem using MATLAB. An 
obvious extension of this work would be to use high-resolution fault grids 
and fast compiled solvers. 
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Figure 25: Inverted Landers slips from the local dataset, only consistent con- 
straints, faults 1-3, /i = 0.7, homogeneous vs layered Earth. 



6 Conclusions and Perspectives 

Results of our study can be summed up as follows: 

• Elastic heterogeneity appears to have a significant qualitative and 
quantitative impact on stress transfer, even when average moduli 
remain constant. 

• Stress transfer is sensitive to the spatial distribution of slips on the 
source faults, and may provide a mechanism for qualitative slip in- 



• Application to Landers and 100 random "aftershocks" , both in a 
neighbourhood of Landers and spread around the entire rupture, 
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Figure 26: Inverted Landers slips from the local dataset, all constraints, faults 
1-3, n = 0.7, homogeneous vs layered Earth. 
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Figure 27: Inverted Landers slips from the decimated (20 faults) regional 
dataset, only consistent constraints, faults 1-3, ^, = 0.7, homogeneous vs layered 
Earth. 
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Figure 28: Inverted Landers slips from the decimated (20 faults) regional 
dataset, all constraints, faults 1-3, /i = 0.7, homogeneous vs layered Earth. 
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Figure 29: Inverted Landers slips using the regional dataset, only consistent 
constraints, faults 1-3, /i — 0.7, homogeneous vs layered Earth. 
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Figure 30: Inverted Landers slips using the regional dataset, all constraints, 
faults 1-3, /i ~ 0.7, homogeneous vs layered Earth. 
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Figure 31: Inverted Landers slips using the Hector Mine event only, un-clamping 
constraints, faults 1-3, /i — 0.8, homogeneous vs layered Earth. 
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Figure 32: Inverted Landers slips from the decimated (20 faults) regional 
dataset, all constraints, high-resolution fault grid, fault 2, /it = 0.7, homoge- 
neous vs layered Earth. 



indicates that, with a few exceptions. Landers appears to have en- 
couraged slips for these events. 

• Application to Hector Mine indicates un-clamping as the dominant 
slip encouraging effect of Landers, in agreement with earlier works. 

• Inverse stress transfer may be a useful technique for qualitative slip 
inversion, however, further study is required to ascertain the extent 
of its applicability and stability to data uncertainty 



7 Appendix A. Reproducibility 

The computational toolkit developed for this study consists of the follow- 
ing components: 

• Elastostatic Green's Tensor Computation Module elastostatic_green.m 

• Surface Displacement Modeling Module model_surf ace_displacements . m 

• Direct and Inverse Stress Transfer Analysis Module transfer _stress . m 

• Auxiliary functions: interp_indices .m, NN_interp_indices .m, ROT.m 
The Green's tensor computation module reads inputs from three files: 

• control parameter file (specified by the $ELASTIC_GT_CTRL environ- 
ment variable in the batch mode or user-specified input); 

• elastic Earth model file ($ELASTIC_GT_MODEL or user input); 

• source faults ($ELASTIC_GT_FAULTS or user input). 

Following are sample input files used for the computation of the Green's 
tensor on the regional Landers/Hector Mine grid for the layered Earth 
model (see above): 

controls .dat : 

nx Sigma ny nz Ix ly Iz 



256 



1.5 



80 



20 



128 



40 



20 



input_hinodel . dat : 



24 



2 21 


0.2 
















5 23 


0.2 
















8 32 


0.2 
















12 36 


0.2 
















15 38 


0.2 
















18 40 


0.2 
















40 45 


0.2 
















Landers . dat : 
















X 


y 


z 


length 


width 


strike 


dip 


rake 


mag 








9 


30 


14 


-5 


90 


180 


6.8035 


22 


-3.9 


9 


25 


14 


-26 


90 


180 


7.1269 


38.5 


-11.7 


9 


30 


14 


-26 


90 


180 


6.9209 



The stress transfer module requires receiver module definition ($ELASTIC_GT_RECEIVERS 
or user input) similar to the source fault file, but excluding the magnitude 
column and including friction and Skempton columns: 



X y 


z 


length 


width 


strike 


dip 


rake 


mu 


B 


50 5.85 


10 


26 


16 


-35 


77 


180 


0.7 


0.5 


44 13.65 


10 


38 


16 


-15 


77 


180 


0.7 


0.5 


16.5 23.4 


10 


26 


16 


-35 


77 


180 


0.7 


0.5 



Constraints of inverse problem ([7| can be limited to consistent ones 
by uncommenting operator continue on line 823 of transfer _stress.m. 
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